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Abstract 

We report on a study of B° —> tt + n~ decays using 29.1 fb _1 of e + e~ an¬ 

nihilation data recorded at the T(45) resonance with the Belle detector at the 
KEKB storage ring. Making no assumptions about the intermediate mechanism, 
the branching fractions for B° —■» D 0 tt + 7t~ and B° —> are determined to 

be (8.0 ± 0.6 ± 1.5) x 10” 4 and (6.2 ± 1.2 ± 1.8) x 10~ 4 respectively. An analysis 
of B° —> candidates yields to the first observation of the color-suppressed 

hadronic decay B° —► D°p° with the branching fraction (2.9 ± 1.0 ± 0.4) x 10~ 4 . 
We measure the ratio of branching fractions B(B° —> D°p°)/B(B° —> D°lu) = 1.6 
± 0.8. 
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1 Introduction 


Exclusive hadronic decay rates provide important tests of models for B me¬ 
son decay [1]. B decays to final states that include a D () or a D*° accom¬ 
panied by two charged pions are interesting, because such decays provide 
a precision testing ground for factorization [2], and because one can search 
for resonant substructure in the final state. At present, only an upper limit 
B(B° —> D°< 1.6 x ICE 3 [3], has been measured. The D ( '*' >0 n + n~ final 
state includes the B° —» ZA*)°p° decay which has not yet been observed [4]. 
As shown in Fig. 1(a) this decay proceeds via an internal spectator diagram, 
and is “color-suppressed” since the color of the quarks produced by the weak 
current must correspond to the color of the c quark and the spectator quark. 
Recent measurements [5] of the branching fractions for the color-suppressed 
decays B° —> D°X°, where A° = i r°, 77 or uj, are all higher than theoretical 
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predictions [6] providing evidence for failure of the naive factorization model 
and suggesting sizable final state interactions (FSI). In the heavy quark limit, 
the QCD factorization model works effectively for color-allowed decays, while 
color-suppressed decays require substantial correction [7]. Assuming SU(3) 
symmetry for the FSI rescattering phase, the observed discrepancy can be 
accommodated and branching fractions, such as B(B° —> D°p°), can be pre¬ 
dicted [8]. It is important to test whether B° —> D°p°, once observed, supports 
the current pattern of QCD effects in color-suppressed B decays. 



Fig. 1. Decay diagrams for B° —> D°p°,D°lu. 

The dominant diagrams for such neutral B meson decays preserve the spec¬ 
tator d-quark and therefore require that the final state neutral light meson be 
produced via its d — d component (Fig. 1(a)). These diagrams predict equal 
branching fractions for D°p° and D°co and for D*°p° and D*°co. Other dia¬ 
grams, such as W-exchange (Fig. 1(b)) or final state interactions can produce 
the u — u state and therefore give different branching fractions. This equality 
is therefore a very sensitive test for small amplitudes in which the spectator 
d-quark does not appear in the final state and the p or u are produced via 
their u — u components [8,9]. 

In this paper, we will report on the study of B° decays that have one D° or 
D*° and two oppositely charged pions in the final state. Inclusion of charge 
conjugate modes is implied throughout this paper. 


2 Data Sample and Selection Criteria 


The data sample used in this paper was collected with the Belle detector at 
KEKB [10]. It is based on an integrated luminosity of 29.1 fb _1 at the T(4S') 
resonance, corresponding to 31.3 million BB events. 

The Belle detector [11] is a large-solid-angle magnetic spectrometer that con¬ 
sists of a three-layer silicon vertex detector, a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cerenkov counters (ACC), a barrel-like 
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arrangement of time-of-flight scintillation counters (TOF), and an electro¬ 
magnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a 
super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect K® mesons 
and to identify muons (KLM). 

Hadronic event selection is described elsewhere [12]. 7 r° candidates are formed 
by combining two photons detected in the ECL, whose invariant mass is within 
a ±16 MeV/c 2 mass window around the 7 r° peak. The 7 r° daughter photons are 
required to have energies greater than 20 MeV. We require the point of closest 
approach to the origin of each track to be within ±5 mm from the beam 
axis and ±3 cm along the beam axis from the interaction point to remove 
background. Tracks identified as electrons (from the responses of the CDC 
and ECL) or muons (from the response of the KLM) are removed. Kaon and 
pion candidates are distinguished by combining the dE/dx information from 
the CDC, time of flight information from the TOF and hit information from 
the ACC. 

D° candidates are reconstructed in the decay modes K~n + , K~7i + n°, and 
Jl ~ 7 r + 7 r“ 7 r + . For D° —> Jl ” 7 r + 7 r°, the 7 T° daughter photons are required to have 
energies greater than 50 MeV and we select regions of the Dalitz plot with 
large decay amplitudes to further suppress the combinatorial background [13]. 
The invariant masses of D° candidates are required to be within 2.5<r of the 
nominal mass. The selected 7 T°s and D °s are then kinematically fit with their 
masses constrained to their nominal values [14]. D*° candidates are formed 
by combining D° and 7 r° candidates and selecting those with mass difference 
8m = M d *o — M d o in the range 0.1400 GeV/c 2 < 8m < 0.1445 GeV/c 2 . 


3 B Meson Reconstruction 


After selecting D° and D*° candidates, we combine them with two oppo¬ 
sitely charged pions to form B candidates. The two oppositely charged can¬ 
didate pions from the B decay are required to come from a single vertex. 
To remove Kg candidates from the sample, the distance of the 7r + 7r _ vertex 
from the beam interaction point in the r — (j) plane is required to be less 
than 0.8 cm. Two kinematic variables are used to identify signal candidates, 
the beam constrained mass, M bc = y (E^ n ) 2 — (PjW) 2 , and the energy dif¬ 
ference A E = Eg M — E£Jf m where E% M and Pg M are the center of mass 
(CM) energy and momentum of the B° candidate, and E™ m = y/s/ 2 = 5.29 
GeV. We select events with |AA| < 0.2 GeV and 5.272 GeV/c 2 < M bc < 
5.288 GeV/c 2 (5.271 GeV/c 2 < M bc < 5.289 GeV/c 2 ) for D°ir + n- ( D*°n+n~) 
final states. Further, if there are multiple B candidates in an event, we choose 
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the candidate with the smallest y 2 combination, 

X 2 = XD° + Xl+ir-i+XSm) (1) 

where, xjjo and y 2 +7r _ are obtained from D° and 7r + 7r _ vertex fitting re¬ 
spectively. For decay modes containing D*°, x 2 m — defined as the square 
of the difference of 8m from its nominal value, in units of its resolution, 
(A(8m)/cr(8m )) 2 — is additionally included in the best candidate selection 
requirement. 


4 Background Suppression 


Since the continuum background (arising from e + e _ —> qq(q = u, d , c, s ) tran¬ 
sitions) has a different event topology, shape variables are very effective at 
improving the signal to noise ratio. Events are required to satisfy i? 2 < 0.5, 
where i? 2 is the ratio of the second Fox-Wolfram moment to the zeroth mo¬ 
ment determined using charged tracks and unmatched neutral showers [15]. 
The angle between the B candidate direction and the thrust axis [16] of the 
rest of the event ( 6t ) is required to satisfy | cos(6U)| < 0.7. 

For the B° —> D W° 7r + 7r branching fraction measurements, we make no as¬ 
sumptions about the intermediate mechanism, except that we reject the large 
contribution from the well-established decay B° —> D* + n~ to the D° 7r + 7r^ fi¬ 
nal state. These events are rejected by requiring M|, 07r+ > 4.62GeV 2 /c 4 (Fig. 
2), which removes 1% of the phase space for B° —» H°7r + 7r^. As the decay 
B° —> H2(2460) + 7t _ is not well established [14], no attempt is made to reject 
it and this mode is thus included in our branching fraction measurement. 

Color-favored decays can also cause a background when a final state pion is 
replaced by a pion from the decay of the other B (for example B~ —> D^°p~ 
may be reconstructed as B° —> ZA*)°7r + 7r - ). To reduce this background we 
veto events which can also be reconstructed in a color-favored mode. This 
requirement removes 1% of the signal candidates. Using a sample of 44 million 
generic b —> c decays generated via Monte Carlo (MC) simulation, the small 
remaining background is studied and found not to peak in M bc or A E. 


5 Branching Fractions for D 0 tc + t: and D*°tt + tt final states 


The distribution in A E for the surviving candidates for B° —» D°h + ti is 
shown in Fig. 3(a). Since intermediate resonances dominate the decay rate we 
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obtain a non-uniform distribution of events on the Dalitz plot. In addition, the 
efficiency varies across the Dalitz plot clue to momentum dependences of the 
reconstruction and particle identification efficiencies. We divide the Dalitz plot 
into six different regions expected to be dominated by different intermediate 
processes as shown in Fig. 2 and determine the efficiency [17] and signal yield 
(from A E fit) for each. Table 1 summarizes our results. 



Fig. 2. Dalitz plot for B° —> D°tt + tt~ events with |AE| < 0.03 GeV, showing the 
regions (a)-(f) used in the efficiency measurement. Events in the dashed box are 
used for the branching fraction measurement of B° —► D°p° as explained in the 
text. 

For each Dalitz plot region we model the signal in A E with a Gaussian func¬ 
tion where both the mean and width are fixed from MG studies. The back¬ 
ground shape in this fit is modeled by two components: (1) a linear shape 
for continuum background obtained from the sideband data (5.20 GeV/c 2 < 
M hc < 5.26 GeV/c 2 ); (2) a smooth histogram shape for B° —> D*°n + 7r~ feed- 
down obtained from MG. The normalizations of the signal and background 








AE (GeV) 

Fig. 3. (a) A E distribution for B —> events satisfying M^ 0 + > 4.62 

GeV 2 /c. (b) A E distribution for B —> D* 0 tt + it~ candidates with no requirement 
on M 2 t(%+ . 

Table 1 

Summary of branching fraction results for B° —r D°ir + ir~ in different regions of the 
Dalitz plot. The last row gives the sums of the signal yields and branching fractions. 


Region Efficiency (%) Signal Yield Branching Fraction (x 10 4 ) 


(a) 

1.87 ± 0.09 

98 ± 15 

1.7 ± 0.3 

(b) 

1.66 ± 0.11 

70 ± 18 

1.3 ± 0.3 

(c) 

1.88 ± 0.08 

17 ± 5 

0.3 ± 0.1 

(d) 

1.94 ± 0.07 

57 ± 15 

0.9 ± 0.2 

(e) 

2.10 ± 0.17 

76 ± 19 

1.2 ± 0.3 

(f) 

1.85 ± 0.12 

150 ± 16 

2.6 ± 0.3 

Total 


469 ± 38 

8.0 ± 0.6 


components are free parameters in the fit. We obtain the branching fraction 
for B° —» Z}°7r + 7r - by taking the sum of the branching fractions in the six 
regions of the Dalitz plot and making a correction of 1% for the unobserved 
region where M |, 0n+ < 4.62 GeV 2 /c 4 . In all branching fraction calculations we 
assume equal production of B°B° and B + B~ pairs from the T(4S'). 

To estimate the branching fraction for B° —» D* q ti + ti~ decays, we make 
no restriction on M^ t0n+ . Due to limited statistics, we do not estimate the 
branching fraction region by region. Instead, we use the yield from the A E fit 
(Fig. 3(b)) and include a model dependent systematic error (19%) that arises 
from the difference between the detection efficiency when the signal MC events 
are B° D*°7c + n~ and B° —> D*°p°. The two detection efficiencies are 0.26% 
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and 0.32%, respectively where the B° —> D*°p° decay is generated with equal 
rates to each hclicity state. 

The background near the lower side of the A E distribution is modeled by 
B + —>• D*°af feed-down measured using MC. The yield from the fit is 62 ± 
12 events. We measure the branching fraction for B° —> H*°7r + 7r - using the 
phase space MC efficiency. The results are summarized in Table 2. 


6 Search for Color-Suppressed B° -> )°p° 

Decays 


Multi-body decays of B mesons can occur through various strong resonances 
that can interfere with each other. We search for color-suppressed B° —» 
£>(*)° p° decays in the Z%*)°7r + 7r final state. We study the 7r + 7r invariant 
mass of the events in the signal region (|Ai?| < 0.030 GeV for D°7r + 7r _ and 
|A E\ < 0.035 GeV for and fit the p° yield with a relativistic Breit- 

Wigner function whose mean and width are fixed to the PDG values [14] to 
estimate the branching fraction. 

To study the color-suppressed decay mode B° —> D°p°, we require M|>o 7r + > 
14.0 GeV 2 /c 4 to remove backgrounds from D* + n~, decays and other 

D resonances. After this requirement, we clearly see an excess at the p° 
mass in the 7r + 7r - invariant mass distribution (Fig. 4(a)). The excess around 
1.45 GeV/c 2 can be modeled by either a p(1450) or an / 0 (1370) resonance; we 
cannot discriminate between these states, or alternative models of the excess, 
based on the fit. Events near 0.5 GeV/c 2 may come from the a [18] resonance. 
We extract the p° yield using a one dimensional likelihood fit. We use a model 
which includes one low mass and one high mass wide resonance. The masses 
and widths are fixed, and the amplitudes and phases are free parameters in 
the fit. The error from the fit therefore incorporates the error from the relative 
phases of the interfering terms: this tends to increase the error on the yield. 
The background under the signal events is described reasonably well by data 
from the A E sideband (0.06 GeV < A E < 0.20 GeV) shown as the hatched 
histogram in Fig. 4(a). We model this shape with a combination of phase 
space, a polynomial and a Breit-Wigner function, where the third term takes 
into account the possible contribution of true p° in the background. From the 
fit, we obtain 86 ± 30 signal events corresponding to a branching fraction 
of B(B° —> D°p °) = (2.9 ± 1.0 ± 0.4)xl0~ 4 . The statistical significance of 
the signal, calculated as ^/—2 ln(Lo/ L max ) where L. max is the likelihood with 
the nominal yield and Lq is the likelihood with the signal constrained to be 
zero, is 6. la. We fold a strong correlation between the amplitude of the p° 
component and its relative phase with respect to the higher mass resonance; 
if the amplitudes and phases of the high and low mass resonances are fixed 
at their obtained values, and the fit is repeated, a p° yield of 86±24 events is 
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obtained. We have repeated the fit with a number of different models includ¬ 
ing vector and scalar resonances at different masses and with different widths; 
the variation in the central value of the p° yield is negligible compared to the 
error from our default fit. As a further cross-check, we examine the the helicity 
angle (0 P ), defined as the angle in the p° rest frame between the direction of 
the 7T + and the p° direction in the B rest frame, and find it to be consistent 
with the expected shape [19]. 

To extract the branching fraction of B° —> D*°p° we require + > 

6.3 GeV 2 /c 4 , which removes backgrounds coming from B° —» Zf2(2460) + 7T - . 
We fit the mass distribution using a relativistic Breit-Wigner function 

after fixing the background shape from the A E sideband (Fig. 4(b)). The 
p° event yield is 29 ± 8, however since the limited statistics prevent us from 
studying possible interferences with other resonances, we cannot interpret this 
as evidence for D*°p° and we set an upper limit of the branching fraction. As¬ 
suming Gaussian statistics, we find B(B° —> D*°p°) < 5.1 x 1CT 4 at the 90% 
confidence level. We summarize our results in Table 2. 



Fig. 4. M n + n - distribution from (a) B° —■> ~ and (b) B° —► D* 0 tt + tt~ final 

states. The histogram represents the data from the A E sideband and the fit to the 
histogram is parameterized as described in the text. 


Table 2 

Summary of branching fraction results for B° —> D ■* )°7T + 7r and B° —> D V [20]. 


Mode 

Efficiency (%) 

Branching Fraction (x 10 4 ) 

Significance ( a ) 

B° — > Z) 0 7r + 7r~ 

1.86 

8.0 ± 0.6 ± 1.5 

18.3 

B° -»• D*°ir+ir- 

0.32 

6.2 ± 1.2 ± 1.8 

6.5 

B° -> D°p° 

0.94 

2.9 ± 1.0 ± 0.4 

6.1 

B° -> D*°p° 

0.24 

< 5.1 

- 
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The following sources of systematic error are considered in our measurements: 
( 1 ) tracking efficiency ( 2 . 0 % per track - measured by comparing the yield of 
the decay modes (i) r/ —> and (ii) rj —> 77 ); ( 2 ) particle identification 

efficiency for n (4.3%); (3) D° reconstruction efficiency and decay branching 
fractions (12.7% - measured by comparing the observed yield of B~ —► D°n~ 
events with the expected yield using known branching fractions [14]); (4) slow 
7T° efficiency (10.7% - measured from the ratio of branching fractions of B~ —» 
D° and B~ —> D*° 7 r _ ); (5) A E signal parameterization (typically 8 %); 
( 6 ) number of BB events (1.0%) and (7) MC statistics (3-5%). As described 
previously, an additional systematic error due to model dependence of the 
efficiency calculation is added for B° —> Z%*)° 7r + 7r - _ The total systematic 
error is obtained by combining the different contributions in quadrature. 


7 Summary 


In summary, we report the first observation of the color-suppressed B° —> 
D°p° decays and measure the branching fraction for B° —> D^°n + n~. Our 
measurement of B(B° —> D°p°) is higher than the factorization prediction of 
0.7 x 10 ” 4 [ 6 ], thus continuing the trend mentioned in the introduction. When 
we compare the branching fraction of B° —>• D°p° to our previous measurement 
of the branching fraction of B° —> D°u [5], we obtain the ratio B(B° —» 
D°p°)/B(B° —> D°u) = 1.6 ± 0.8. The error includes both statistical and 
systematic errors where the correlation of the systematic errors has been taken 
into account. Future measurements with more statistics will allow precise tests 
of the mechanisms involved in color-suppressed B decays. 
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